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Abstract. We give the Fuchsian linear differential equation satisfied by x' 4 '> 
the "four-particle" contribution to the susceptibility of the isotropic square lattice 
Ising model. This Fuchsian differential equation is deduced from a series expansion 
method introduced in two previous papers and is applied with some symmetries 
and tricks specific to x' 4 ' ■ The corresponding order ten linear differential operator 
exhibits a large set of factorization properties. Among these factorizations one is 
highly remarkable: it corresponds to the fact that the two-particle contribution 
X^ 2 ' is actually a solution of this order ten linear differential operator. This result, 
together with a similar one for the order seven differential operator corresponding 
to the three-particle contribution, x i leads us to a conjecture on the structure 
of all the n-particle contributions x ■ 
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1. Introduction 

The magnetic susceptibility of square lattice Ising model, can be written [1] as an 
infinite sum 



of individual contributions, with the odd (respectively even) n corresponding to high 
(respectively low) temperature case. These individual contributions are [n — 1)- 
dimensional integrals [2, 3, 4, 5, 6, 7], and are seen as successive n— particle 
contributions to the susceptibility [1]. 
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To get an understanding of the analytical structure of \, two approaches are 
usually taken nowadays. One approch taking into account a fundamental non-linear 
symmetry, namely non-linear Painleve difference equations [8, 9, 10, 11, 12, 13]. 
provides a series expansion for the whole susceptibility x- With this method 
coefficients of x were recently generated [8] . 

The second approach considers the individual n-particle excitations as given by 
(n— l)-dimensional integrals. Isotropic series coefficients are generated [6, 7, 8]. This 
latter method allows one to seek the differential equations satisfied by the x^ n \ since 
they are D-finite contrary to the whole susceptibility x f° r which there are strong 
indications that it has a natural boundary in the complex plane of the variable 
s = sh(2K), where K = J/kT is the usual Ising model coupling constant. This 
has been shown for the isotropic case [6] and for the anisotropic case [14, 15, 16]. Such 
a function cannot be Infinite. 

The understanding of the magnetic susceptibility may then require the knowledge 
of each (or some) of the individual contributions. This knowledge can be in the form 
of a closed expression, as is the case for and x^ or in the form of a differential 
equation as found for x^ [17, 18]. The last case was far from being obvious and has 
required the building of an original method of expansion. The use of a remarkable 
formula allowed us to give the series expansion in the temperature variable (or a closely 
related variable) where the (n — l)-dimensional integrals have been fully performed. 

In this paper, we continue to use this expansion method to tackle the next 
individual contribution, namely, x' 4 '- We should note that, although, the method 
is general and applicable for high or low temperature and for any n, some of the tricks 
and tools used may be specific for a given x^ ■ I n Section 2, we present the basic 
features of the expansion method that allow us to obtain the fully integrated x^ 
as four sums of products of four hypergeometric functions, without any numerical 
approximation. In Section 3, we give the homogeneous Fuchsian linear differential 
equation satisfied by x' 4 ''- Section 4 contains some remarkable algebraic properties of 
this differential equation. Finally, Section 5 contains our conclusions. 

2. Fully integrated x' 4 ' expansion 

2.1. The expansion method 

Let us focus on the fourth contribution to the susceptibility % defined by the triple 
integral as given in [7] 
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Instead of the variable s, we found it more suitable to use w — ^s/(l + s 2 ) which 
has, by construction, Kramers- Wannier duality invariance (s «-> 1/s) and thus allows 
us to deal with both limits (high and low temperature, small and large s) on an equal 
footing [17, 18]. It is also convenient to consider the scaled variables 

2 (8) 
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which behave like 1 + 0(w) at small w. 

As performing the integrals in (2) is highly non trivial, we apply the "expansion 
method" previously described in [17, 18], where the key ingredient was the Fourier 
expansion of yx n , a quantity appearing in any ■ 

This remarkable formula, for yx n , that carries only one summation index reads 



yx n = a(0,n) + 2 w k a(k, n) cos k<j> 
fe=i 

oo oo 
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with 



A(k,n) = A(-k,n) = w w a(\k\ , n) (10) 
where a(k,n) is a non terminating hypergeometric series that reads 

a(k,n)=(^y (11) 

„ /(l + m) (1 + m) (2 + m) (2 + m) _ , , 9 \ 

4^3 ( l 2 - , V 2 % V 2 7 , V 2 - ; 1 + fc, 1 + n, 1 + to; WJ 

where to = fc + n. Note that a(fc, n) = a(n, fc). 

The integrand of x' 4 ^ is expanded in the various variables Xj, instead of the 
variable w. In this framework, with the help of the Fourier expansion (9), the angular 
integration becomes straightforward as was shown in [18] for the case (see below 
for the x' 4 ^ case). 



2.2. Calculation o/x (4) 

With taken as m (4), x (4) will be given by five summations on products of four 
hypergeometric functions. This is shown in Appendix A. This route is feasible for any 
X*-™-* and does not use any symmetry or tricks specific to the considered x'™'- 

In the following, we use alternatively a simplified form of equivalent, for 
integration purposes, to (4) such that x' 4 ^ will be expressed by only four summations 
(i.e., one summation less compared to the expression (88) given in Appendix A). Let 
us just sketch the salient steps of this calculation. The details are left to Appendices 
B, C and D. 
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the quantity becomesf : 
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Using the symmetry of the rest of the integrand in the angular variables, the quantity 
can be written (see Appendix B) asj 

-!<*-*><* -^<™?- 

(Pft-s. - (1 « 3) - (2 « 3)) (14) 

where : 

P 12-34 = (Zl - X 2 )(X 3 - X A ) ■ (15) 

/ x 1 x 2 x 3 x 4 (x!X 2 ) 2 {x 3 x 4 ) 2 \ 

V 1 — ^1^2 1 — ^3^4 1 — ^1^2 \ — X?,XA> 

Taking the expression (14) for H^, expanding the integrand in the Xj variables, one 
notes that the integrand depends only on combinations of the form 

(f[yi ■ x?J ■ (Zi Z 2 )(Z 3 Z A ) (1 Zl ^ f (16) 

which have simple integration rules (see Appendix C) and, thus, the problem of angular 
integrations is settled. 

Finally, the expansion method described above, together with the form (14) of 
H( 4 \ allow us to obtain X^ 4 \ w ) as a fully integrated expansion in the form (see 
Appendix D) 
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2 

where 9(x) is the step function defined as 



«- r) =<l: ill w 



and 



V(m, k, n,j) = d(m, m + l;j + k + 2) • d(n + m + 2, n + m + 3; k) 
+ d(m, m + 1; k) ■ d(n + m + 2, n + m + 3; j + k + 2) 
+ d(m, n + m + 3; j + k + 2) • d(m + 1, n + m + 2; k) 
+ d(m, n + m + 3;k) ■ d(m + 1, n + m + 2; j + k + 2) 

— d(m, n + m + 2;j + k + 2) • d(m + l,n + m + 3;k) 

- d(m, n + m + 2; k) ■ d(m + 1, n + m + 3;j + k + 2) (19) 

t We have used the constraint (6) in the form FJ^ Zi = 1. 

J Throughout this paper the notation = stands for equality for integration purposes. 
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with : 

d(ni,n 2 ; k) — a(ni, k + 1) a(n 2 , k) — a{m,k) a(n 2 ,k + 1). (20) 
2.3. Series generation 

Note that the summand in equation (17) depends only on combinations of 
Cd(ni,ri2,k,j) which is the coefficient of w 2 ^ in the expansion of d{n\, n 2 ; k). It 
is given by 

i 



Cd(ni,n 2 ,k,j) = ^2 (c a (ni,k + l,i)C a (n 2 ,k,j - i) 

i=0 

-C a ( ni ,k,i)C a (n 2 ,k + l,j-i)} (21) 
where C a (n, k, i) is the coefficient of w 2t in the expansion of a(n, k) which reads : 

c„(n,M>=(" + r ')(":::*) <*> 

From our integrated form of x' 4 ', the generation of series coefficients becomes 
straightforward. Recall that (17) is already integrated, and, thus, the computing time 
to obtain the scries coefficients comes from the evaluation of the sums. For the forms 
used, this time is of order N 7 . Improvements can be made by optimal data storage 
in order to avoid repeated summation evaluations. Actually, we have found it more 
efficient to store the coefficients C<j. 

We have been able to generate, from formal calculations, a long series of 
coefficients from the expression (17) up to order§ 432 : 

* 6 ^ 16 = 1 + 64 W 2 + 2470w 4 + 74724 w 6 + ... + 0{w 434 ) (23) 

Note that our formal calculation program has been rewritten by J. Dcthridgc|| into 
an optimized C++ program that can give the series in few hours using very little 
memory. 

3. The Fuchsian differential equation satisfied by \^ 

It is clear from the expression (17) that x^is even§ in w. We thus introduce, in the 
following, the variable x — 16 w 2 . With our long series, and with a dedicated program, 
we have succeeded in obtaining the differential equation for x' 4 ^ that is given by (with 
x = 16 w 2 ) 

10 d n 

E a ^ x ) ■ d^ F ^ - ( 24 ) 

§ More precisely we obtained this series up to order 368 with a 2 Giga-memory computer, up to 
order 390 with a 3 Giga-memory computer, and up to order 432 with a 8 Giga-memory computer of 
the stix laboratory at the Ecolc Polytechnique (mcdicis platform for formal calculations). 
|| Private communication. 

§ This is also the case for any x' 2 ™' (see for instance [7]). 
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with 

aio = -512 x e {x - 4) (1 - xf P w (x), 

a 9 = 256 (1 - xf x 5 P 9 (x), a s = -384 (1 - xf x 4 P 8 (x), 

a 7 = 192 (1 - xf x 3 P 7 (x), a 6 = -96 (1 - xf x 2 P 6 (x), 

a 5 = 144 (1-x) xP 5 (x), a 4 = -72P 4 (x), a 3 = -108P 3 {x) 

a 2 = -54P 2 (x) ai = -27Pi(x) (25) 

where Pio(x), P$(x) ■ ■ ■, Pi(x) are polynomials of degree respectively 17, 19, 20, 21, 
22, 23, 24, 23, 22 and 21 given in Appendix E. 

With (25), the differential equation needs 242 unknowns to be found (counting 
the polynomial in front of the derivative of order 0, which is identically null). Our 
series expansion for x^/16w le , having only 217 terms in the variable x = I6w 2 , is, 
thus, not long enough to let the differential equation be found. This calls for some 
comments. 

The differential equation (24) is of minimal order. It is obvious that one can 
obtain other differential equations of greater order. As explained in [18] (see Section 4), 
before the differential equation built from a series expansion pops out, the singularities 
computed as roots of the polynomial in front of the highest derivative reach stabilized 
numerical values as the degrees of the polynomials and/or the order get higher. Even 
the minimal multiplicity of the singularity can be seen. This leads us to take the 
following form for the differential equation 

Ri(x) ■ x* (1 - xf ■ £rF(x) = (26) 

where only the physical singularities are explicitly included. Note that the other (non- 
apparent) singularity can also be included. If the differential equation can be identified 
with the number of series coefficients at hand, the unnecessary terms will factor out. 
A compromise has to be found with respect to how long the series is. Now, there are 
many ways to choose the degrees of the polynomials Ri{x). We have taken 

deg(Ri{xfj = v + q-i, i = Q,l,---,q (27) 

in order to have the point at infinity as a regular singular point\\. 

If the differential equation of order q and degree fi exists, the series of x (4) /!6w 16 
should have at least N terms, with : 

N = M .(g+l) + i g .(g + 3 ). (28) 

Staying below the above hyperbola, we have obtained, at q = 11 and /i = 9 which 
requires 185 terms in the series, two linearly independent differential equations that 
satisfy the remaining 30 terms of the series. The combination of these two differential 
equations gives (24) which has been checked to be of minimal order. Note the fact 
that besides the Fuchsian differential equation of minimal order with an apparent 
polynomial, there are other differential equations of higher order that require less terms 
in the series to be identified. We plan to report on this feature elsewhere. 

The singularities of this differential equation x = 0, 1, 4 and x — oo are all regular 
singular points. The roots of the polynomial P\o{x) are apparent singularities. One 

|| This feature can easily be seen, making the change of variable t = l/x and looking for the necessary 
condition for t = to be a regular singular point from which (27) is deduced. 
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notes, contrary, to x^ that no new singularity is found besides the known physical 
and non physical singularities (i.e., Nickel's [6, 7]). The critical exponents of all these 
singular points are given in Table 1 with the maximum power of logarithmic terms 
in the solutions. As was the case for x^ , these logarithmic terms appear due to the 
multiple roots of the indicial equation. A noteworthy remark is the occurrence of 
logarithmic terms up to the power 3 for x^ to be compared with the power 2 for x^ 
at the singular points w = ±1/4 and w = oo. 



x- singularity 


s-singularity 


Critical exponents in x 


P 





0, oo 


8,3,3,2,2,1,1,0,0,-1/2 


1 


1 


±1 


3,2,1,1,0,0,0,0,-1,-3/2 


3 


4 




8,7,13/2,6,5,4,3,2,1,0 





oo 


±i 


5/2,3/2, 3/2, 1/2, 1/2, 1/2, 1/2,0, -1/2, -1/2 


3 


xp, 17 roots 


sp, 68 roots 


10,8,7,6,5,4,3,2,1,0 






Table 1: Critical exponents for each regular singular point. P is the maximum 
power of the logarithmic terms for each singularity, xp is any of the 17 roots of Pio(x). 



It is worth recalling the Fuchsian relation on Fuchsian type equations. Denoting 
by x\, X2, x m , x m+ i = oo, the regular singular points of a Fuchsian type 
equation of order q and Pj,i, Pj.q (j = 1, ••• ,m+ 1) the q roots of the 

indicial equation [28, 29] corresponding to each regular singular point Wj, the following 
Fuchsian relation [28, 29] holds : 

m+i q , i\ / i\ 

E E - (m ~ 1} 2 g( ■ (29) 

j = l fe = l 

The number of regular singular points here is m + 1 = 21 corresponding respectively 
to the 17 roots of Pio, the x = 0, 1, 4, regular singular points, and the point at 
infinity x = oo. The Fuchsian relation (29) is actually satisfied here with q =10, 
to = 20. 

Considering the Fuchsian relation (29), but now in the variable s, one first 
remarks that the p 3i k exponents are the same as the ones in the x variable for the 68 
roots of the apparent polynomial, as well as for the four x = 4 singularities (namely 
s = ±5 ± i^f-), but are multiplied by a factor two for all the other ones. This is 
obvious from the definition of x = 4s 2 / (1 + s 2 ) 2 at x — and x = oo, and this can be 
seen for x = 1 from : 

(30) 

The Fuchsian relation is actually satisfied in the s variable, but now, with q = 10, 
to = 77. 
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Let us close this section by confirming the dominant singular behavior at x = 4 
and x = oo given in [7] , since this can easily be done from the critical exponents given 
in Table 1. Near x — 4, the solution behaves as i 13 / 2 and near x = oo as i" 1 / 2 log(i), 
with respectively t = 4 — x and t — 1/x. Fom Table 1, the subdominant singular 
behavior at x = oo is i 3 / 2 log(i) and t 1 / 2 \og k {t) (with k = 1, 2, 3). 



4. Properties of the Fuchsian differential equation (24) 

From Table 1, and from the formal solutions around the singular points, it is easy to 
find the following simple solutions of (24) : 

So(x) = constant (31) 

„ , . 8-12x + 3x 2 

8(1 -,)3/ 2 * ( 32 ) 

52W= 2(1-,)^ - (33) 
These solutions correspond to solutions of some differential operators of order one. 
Let us call these order one differential operators respectively L , L\ and L 2 . 

A remarkable finding is the following solution of the Fuchsian differential equation 
(24) : 

Ss{x)= h x2 ' 2Fi (l'l ]3]X ) (34) 

which is nothing but the two-particle contribution to the magnetic susceptibility, i.e., 
associated with an operator of order two, N (iVo(<5>3) = 0). We will come back 
to this point later. 

The second solution of the order-two operator No is given in terms of the MeijerG 
function [19] : 

Ss(x) = | MeijerG ([[], [1/2, 3/2]], [[2, 0], []},x) (35) 

which can also be written as 

S 3 (x) = S 3 (x) log(x)+ B(x) (36) 

with : 

With these five solutions corresponding to three differential operators of order 
one, and one differential operator of order two, it is easy to construct 24 factorizations 
of £io, the differential operator corresponding to the Fuchsian differential equation 
(24), which can be written^ as : 

C w = 5 ■ G(N). (38) 

G(N) is a shorthand notation of a differential operator of order 5, factorizable in 
one operator of order two and three operators of order 1. G(N) has 24 different 
factorizations involving eight differential operators of order two and 24 operators of 

K All the operators are such that the coefficient in front of the highest derivative is +1. 
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order one. The differential operator 5 factorizes as M x -£24, i.e., one operator of order 
four and one operator of order one. Let us give two examples of the 24 factorizations" 1 ": 

£10 = Mi • £ 24 ■ N A ■ L 12 ■ £3 • £ (39) 
£10 = Mi • £ 24 ■ £13 -N 6 -L 3 -L a (40) 

This large number of factorizations [22, 31, 21] induces the occurrence of 
intertwiners f. In the examples above, one has N4 ■ L\ 2 — £13 • Nq. Seeking a 
similar relation for £24 ■ N4, one finds Ng ■ £25. This last factorization introduces six 
factorizations that we denote as: 

C w = Mi • N 9 ■ G(L). (41) 

G(L) is a notation for an operator of order four, factorizable as four operators of order 
one. 

One factorization of G(L) reads: 

G(L) = £ 25 • £12 • £3 • £o- (42) 

This differential operator G(L) that factorizes £10 at right, obviously has So, Si and 
S 2 as solutions. The fourth solution (of the order four differential operator G(L)) can 
be obtained by order reduction. It reads 



4(z-2)V4— ^ 
04(2;) = h 16 log 



(2 + V4 - x) 2 

+ WS 1 (x)-\ogg(x) -16V*~S 2 (x) . ^ (I, I; 1 1) (43) 



with : 



g(x) = i ((8 - 9x + 2x 2 ) + 2 (2 - x) ■ - x){A - xf) (44) 

To get more factorizations, we need to obtain simple solutions of £^ , the adjoint 
of the differential operator £10 • There is no solution of C* w corresponding to an 
order-one operator, however, we have been able to find a solution corresponding to an 
operator of order two (denoted Ng). This solution! of Ng is a combination of elliptic 
integrals with polynomials of quite large degrees and reads 

K(x)= 2 F 1 (±,±;l;x), E(x) = 2 F Y (-1 1; 1; x) (46) 

qx (x) = (1 - a;) (47352014438400- 246257318625280 x 
+ 275880211382272 x 2 + 68328139784192 x 3 
+ 645943284072448 x 4 - 2774821715853312 x 5 
+ 3217221650489344 x 6 - 683914539437568 x 7 

+ We denote by L the operators of order 1, by TV, the operators of order 2 and by M the operators 
of order 4 (see Appendix F). 

f We thank Jacques-Arthur Weil for usefull comments on the equivalence of linear differential 
operators [20] (see equation (5) in [21]). 

X We did not look for the second solution of Ng, our purpose being the factorization of £10. 



where 



and : 
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- 2042467767948624 a; 8 + 3083863919521506 a; 9 

- 2746206480894969 a; 10 + 1558224994851490 a; 11 

- 347705392468761a; 12 - 145625559012638 a; 13 
+ 117842186745065 a; 14 - 30744722745590 a; 15 
+ 3089482306025 a; 16 + 18651488480 a; 17 

- 21574317760 a; 18 + 821760000 a; 19 ) (47) 

q 2 ( x ) = - 47352014438400+ 269933325844480 a; 

- 390130367987712 a; 2 + 27877970018304 a; 3 

- 580571855978496 a; 4 + 3135069528473600 a; 5 

- 4488375407386624 a; 6 + 1736922901371392 a; 7 
+ 2600277912748368 a; 8 - 5123144341863018 a; 9 
+ 5224617790090830 a; 10 - 3547418998359865 a; 11 
+ 1453586336314895 a; 12 - 273126255420088 a; 13 

- 8194519962996 a; 14 + 11308926014655 a; 15 

- 1235672485785 a; 16 - 60982101700 a; 17 
+ 15765744320 a; 18 - 607865600 a; 19 

- 3840000 a; 20 (48) 



This order 2 operator N$ completes the factorization scheme of Ciq, the 
differential operator of order ten. Let us recap the factorizations 



G(L) is a differential operator of order four, factorizable in order one operators. 
It has six different factorizations involving thirteen differential operators of order one. 
G(N) is an operator of order five, factorizable in one operator of order two and three 
operators of order one. G(N) has 24 different factorizations involving eight differential 
operators of order two, and 24 order one operators with twelve appearing in G{L). 

All these 36 factorizations (49) are given in Appendix F. Appendix F shows that 
these 36 factorizations can be considered as only one factorization up to a set of 
equivalence-symmetries [20, 21, 22, 23, 24, 25, 26]. 

From the ten solutions of the differential equation (24), six solutions are given 
explicity, So, Si, S2, S4 (this last being a solution of an order four operator) and two 
solutions (£3, 1S3) corresponding to the differential equation of \^ ■ The remaining 
four solutions are those of the operator of order eight, M 2 • G(L). 

From the 36 factorizations shown in Appendix F, those six of the form N% ■ M 2 ■ 
G(L) are of the most importance. Their occurrence allows us to get the contribution a 
of %( 2 ) in the physical solution \^ from M 2 -G(L)(x^ — a\^) = 0. This contribution 
is obtained easily and gives 



where $4 is solution of the order eight differential operator M 2 ■ G{L). 



£10 - A^s • M 2 ■ G(L) 
£10 = Mi ■ N 9 ■ G(L) 
£10 = Mi ■ L 2A ■ G(N) 



(49) 



X (4) = ^X (2) + 4>4 



(50) 
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Recall that the same situation occurred for the differential equation of order 
seven for . One can see that the rational solution w/(l — Aw) occurring§ in the 
differential equation of order seven for is nothing but ■ With this remark we 
can rewrite a decomposition of x*- 3 ^ we gave in [17, 18] as follows : 

X (3) = \x {1) + <J>3- (51) 
where $3 is solution of the order six differential operator (noted L 6 in [17, 18]). 



5. Comments and speculations 

Denoting L4 the order ten differential operator associated with the ordinary 
differential equation satisfied by x^ 4 -* > and more generally, L„ the differential operators 
associated with the ordinary differential equation satisfied by x^ , one has : 

L 4 (x (4) ) = U(X {2) ) = 0. (52) 

Furthermore, one can see that the rational solution w/(l — 4w), occurring in the 
differential equation of order seven for x^\ is nothing but y^ and thus one has : 

L 3 (x (3) ) = L 3 (x (1) ) - 0. (53) 

Both relations come from the factorizations of the differential operators corresponding 
to x (3) and x (4) - 

At this point, it is tempting to make some conjectures generalizing relations (52, 
53) and relations (50,51). 

From (52) and (53) the conjecture is 

L 2 „+i(x (2n+1) ) = Wi(x (1) ) = 0, 

L 2 „(x (2n) ) = Wx (2) ) = (54) 

meaning that the differential operator of x^ (respectively x^) right divides the 
differential operator corresponding to x^ 2n+1 ^ (respectively x^ 2 ™')- 

One stronger conjecture is to expect the same situation as in (50,51) occurring in 
the higher particle contributions, i.e., 

= a 2 „x (2) +$2n, 

X (2 " +1) = «2„ + lX W + $2n + l (55) 

where $ 2 „ (respectively $ 2n +i) is solution of a differential operator that right divides 
the differential operator L 2 „ (respectively L 2 „ + i) and is not divisible by the differential 
operator L 2 (respectively Li). In this situation, it is easy to obtain the numbers a n 
which give the contribution of x^ and x' 2 ) in the higher x^'s as explained for x^ 
(see text before (50)). 

Let us note that both conjectures are free from any constraint due to the 
singularities that occur in the differential equations, since any x( 2 ™ +1 ) (respectively 
X (2n) ) has the singularities occurring in x^ 1 ) (respectively x' 2 ')- 

A much stronger conjecture is to expect any x^ to be "embedded" in any x^ 
where m divides n with same odd-even parity. 

Mx (n) ) = Mx (m) ) = 0, 

L„ = 4 m ) -L m . (56) 

§ This corresponds to the solution Si given in [17, 18] which is Si = x (1) /2- 
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This means that L n might be built from the least common left multiple (lclm) of 
the differential operators associated with the L m 's where m divides the integer n 
respecting even-odd parity. The fact that n should be multiple of m is due to the 
non physical singularities appearing in the differential equations. For instance, no 
"embedding" like being solution of the differential equations of X*- 5 -* has to be 
expected. Whether x^ is a solution of the differential equation satisfied by, e.g., 

is not ruled out, the (non-apparent) singularities of the first being (non-apparent) 
singularities of the last. But, then, the new singularities discovered [17] for x^ have 
to occur in the differential equation of x^ ■ 

Similarly one can expect, with this last conjecture, the following relations for Li 2 
and Li 5 (associated with x^ 12 ) and x^) 

L 12 ( X (12) ) = L 12 (x (6) ) = L 12 ( X (4) ) = L 12 ( X (2) ) = 

L 15 (x (15) ) = L 15 ( X (5) ) = L 15 ( X (3) ) = L 15 ( X (1) ) = (57) 

and thus : 

T t (6) t _ T (4) T t (2) -r 

L15 = Li5^ • L5 = ■ L3 = L^' - Li. (58) 

For instance, finding the differential operator satisfied by x*- 1 ^ would enable us to 
see if the differential operator L15 actually has some relation with the least common 
left multiple of Li, L 3 and L 5 . 

6. Conclusion 

Considering the isotropic Ising square lattice model susceptibility, we extended our 
"expansion method" (that allowed us to find the differential equation satisfied by x^ 3 -*), 
to the four-particle contribution to the susceptibility, namely x^ ■ We first obtained 
X^ as a fully integrated multisum on products of four hypergeometric functions, and 
obtained a long series for x^ • From this long series, we gave the Fuchsian differential 
equation of order ten satisfied by x' 4 ^ • This differential equation has a rich structure 
in terms of factorizations. 

We have given in closed form six, of the ten solutions of the differential equation 
(24), namely, So, Si, S 2 , S4, S3 and S3. The remaining four solutions are those of an 
operator of order eight. 

One of these solutions is highly remarkable: it is actually the two-particle 
contribution x^ ■ A similar situation also occured for the differential equation of 
the three particle contribution x^ 3 \ which actually had x^ as solution. 

In general, for all the x^' s : it is tempting to expect x^ to be a solution of 
the differential equation satisfied by the x' 2n+1 - ) ' s and x^ to be a solution of the 
differential equation satisfied by the x^ 2 "^'s. 

Beyond, one can contemplate much stronger conjectures corresponding to further 
"embedding", namely x ( '™' ) being solution of the differential equation of x'™\ when n 
is a multiple of m, n and m having the same parity. 

A confirmation of these various embeddings and conjectures is crucial because it 
corresponds to a global structure of the X^' s an d thus of x itself. We will investigate 
this feature in a future publication. 
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7. Appendix A 

In this Appendix, we show that, with taken as in (4), i.e., not the alternative 
simplified form we use in this paper, our expansion method applied to \^ produces 
five summations on products of four hypergeometric functions. Note that this method 
is applicable to any x^ ■ 

Let us write the product ■ as 

fl(4) . H W = T (4) . A {A) (5Q) 

with : 

T (4) = R (4) , TT {m) 
(1 - x t x J ) 2 

A{i) = i ■ - ^) 2 - ( 6i ) 

i<j 

By standard expansion of the quantity defined by (60), in the variables Xi, 
one obtains 

oooooooooooooo 6 

r (4) = EEEEEEE( 1 +%- 1 )MI^ + 1 )- 

p—0 ii—0 22—0 is—O 1 4—0 «5— Iq— k—1 

~p+ii+22+*3+3 ~p+i2+U+*6+3 ~p+i3+*5+«6+3 

X-^ Jj<2 Xg X^ \\jZiJ 

where 9{x) is the step function given in (18). Defining 

m = p + ii + i 2 + is, n 2 = p + ii+ii + ih 

n 3 = p + i 2 + i4 + ie, "4 = P + h + h + (63) 
one can solve (63) in the indices (is, 14, i$, i%), to obtain: 

is = - (p + h +12) +ni (64) 

U= - (p + ii + i 2 ) + - (ni + n 2 + n 3 - n 4 ) (65) 
»5 = «2 - ^ («i - n 2 + n 3 - n 4 ) (66) 

k = k - ^ ("1 + n 2 - n 3 - n 4 ) (67) 

All these indices should be integers, inducing constraints on the n^'s and limitations 
on the summations of the remaining indices (i.e., p, ii, i 2 ). With this change of 
summation indices, becomes 



00 00 00 00 



E E E E 

7ii—0 ri2— ns— n^—O 

C(n u n 2 ,n 3 ,n A ) • i" 1+3 ^ 2+3 i™ 3+3 ^ 4+3 (68) 
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and, since is completely symmetric in the variables x i7 the coefficient 

C{n\,n 2 , n 3 , n 4 ) is completely symmetric in the n; indices and is given by summation 
on the remaining indices p, i\ and i 2 : 

GO GO go 6 

C(n u n 2 ,n 3 ,n 4 ) = EE E (1 + 6(p - 1)) • J](i fe + 1) 

p=0ii=0i 2 =0 fe=l 

x 0(i 3 ) • 0(i 4 ) • 0(i 5 ) • 0(i 6 ) • <r(ni + n 2 + n 3 + n 4 ) 



'(^ (ni + n 2 + n 3 - n 4 )) (69) 



where theta functions (18) of indices (i 3 , i 4 , i^,io) take place to keep track of the fact 
that these indices should be positive integers. The symbol cr(n) defined as: 

*(n) = ±(l + (-l) n ) (70) 

comes from the fact that i 3 , i 4l 25, ie are integers, due to the right most terms at right- 
hand-side of (65-67), the n^'s should verify n x + n 2 + n 3 + n 4 = even integer. The 
argument in the last theta in (69) comes from (65). 

Let us define the index q and its upper limit of summation (see (64,65)): 

q =p + h+i 2 (71) 

ni +n 2 +n 3 - n 4 
q = mm(ni, - ) (72) 

Furthermore, the coefficient C(ni, n 2 , n 3 , n 4 ) being symmetric in all its 
arguments, it is sufficient to compute it in the case where: 

n\ < n 2 < n 3 < n 4 (73) 

The coefficient C(ni, n 2 , n 3 , n 4 ) now becomes: 

C(ni, n 2 , n 3 , n 4 ) = 0{ ) • <r(ni + n 2 + n 3 + n 4 ) 

qo q / q-p 6 \ 

EE (i+%-i))-E(n^ + 1 )) ( 74 ) 

5=0 p=0 V ii=0 k=l / 

where the indices ? 3 ,i 4 , 25,^6 are given by (64-67) and i 2 = q — p — i\. The three 
summations can be performed and one obtains 

C{n u n 2 , n 3 , n 4 ) = 6{ - ) • cr(m + n 2 + n 3 + n 4 ) 

'qo + 4 

with the polynomial Q(ni, n 2 , n 3 , n 4 , go) given by 

Q(n 1 ,n 2 ,n 3 ,n 4 ,q ) = ^Qo(ni,n 2 ,n 3 ,n 4 ) + -^^Q 1 (n 1 ,n 2 ,n 3 ,n 4 ) 

+ Y^Q2{ni,n 2 ,n 3 ,n 4 ) + Y^Q 3 (n 1 ,n 2 ,n 3 ,n 4 ) 

+ ^ ( 92 - 81n i - 9n 2 - 9n 3 + 63n 4 ) + ^| (76) 
541) 135 



Q{ni,n 2 ,n 3 ,n 4 ,q ) (75) 



with 



Qo(ni,n 2l n 3 ,n 4 ) = (1 + ni)(ni + n 2 - n 3 - n 4 - 2) 

(ni - n 2 + ra 3 - n 4 - 2) (2 + n x + n 2 + n 3 - n 4 ) (77) 
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Qi{n 1 ,n 2 ,n 3 ,n 4 ,) = -2216 + 2472m + 1194n^ - 1281n? + \2&n\ 

- 48n 2 + 1392nin 2 - 987nfn 2 + 126n?n 2 + 774^ + Ulnml 

- \2<on\n\ + \AHn\ - 126nm^ - 48n 3 + 1392nin 3 - 987n?n 3 
+ 126njn 3 - 1548n 2 n 3 - 882nin 2 n 3 + 252nin 2 n 3 - 147n 2 n 3 
+ 126nin2n 3 + 774ng + AAln^l - \2§n\n\ - 147n 2 n 3 

+ \2§n x n 2 n\ + \Mn\ - Y2&n x n\ - 2568n 4 - 708nin 4 

+ 2709njn 4 - 378njn 4 - 288n 2 n 4 + 1134nin 2 n 4 - 252njn 2 n 4 

- 147n 2 n 4 + 126mn 2 n 4 — 288n 3 n 4 + 1134n 4 n 3 n 4 — 147n 3 n 4 
+ 294n 2 n 3 n 4 — 252nin 2 n 3 n 4 — 2b2n\n 3 n 4 + 126n\n 2 n\ 

+ \2§n 4 n\n 4 - 486?t,4 - \hlhn 4 n\ + 378n^n 4 - 147n 2 n 4 

- lA7n 5 n\ + 126mn 3 nl + lA7n\ - \2§n Y n\ (78) 

Q 2 {n 1 ,n 2 , n 3 , n 4 ) = -2992 - 4446n 4 + 6534n? - 1449n? - 1926n 2 

- 819n?n 2 - 270^ + 945nm^ + 315^ - 1926n 3 + 2808nin 3 

- 819nfn 3 + 540n 2 n 3 - 1890nin 2 n 3 - 315njn 3 + 2808n!n 2 

- 270^ + 945nin^ - Zlbn 2 n\ + 315n^ + 594n 4 - 315n 3 n^ 

- 8532n 4 n 4 + 3213nin 4 - 1728n 2 n 4 + 1134nin 2 n 4 + 315n| 

- 315n 2 n 4 — 1728n 3 n 4 + 1134n 4 n 3 n 4 + 630n 2 n 3 n 4 

- 315n^n 4 + 1998n^ - 2079nm 2 4 - i\hn 2 n\ (79) 

Q 3 (ni, n 2 , n 3 , n 4 ) = 520 - 1215n 4 + 495n^ - 207n 2 + 144nin 2 
+ 144nm 3 + 270n 2 n 3 - 135n 3 + 801n 4 - 738nm 4 

- 108n 3 n 4 + 243^ - 135^ - 207n 3 - 108n 2 n 4 (80) 
From the definition of x^ in (2), and from (59,68), one gets : 

oo oo oo oo 

X (4) = 16 ™ 16 • w ni+n2+n3+ni X (81) 

rii— U2— Ti3— ri4— 

C(ni,n 2 ,n 3 ,n 4 ) • M(ni, n 2 , n 3 , n 4 ) 

The expansion method on the variables Xi introduces some summations free of any 
angular dependence. These summations appear in the coefficient C(n\, n 2 , n 3 , n 4 ) 
which is a kind of "geometrical factor" that appears for all the X^' s with n > 4. 

At this step the integrations have not yet been performed. They are contained 

in: 

M(ni ' n2 ' n3 ' n4H 7„ Wo 25T 

/ 4 \ 



Y[ yi -x? +3 )-AW (82 ) 



=i 



Due to its structure, it is obvious that M(ni, n 2 , n 3 , n 4 ) is completely symmetric 
in the indices n^. We write M(ni,n 2 ,n 3 ,n 4 ) as an integral over the four angles 
(0i, 02, 03, 4 ), by introducing Dirac delta function S((f>i + 2 + 3 + 4 ) as a Fourier 
expansion that reads: 



oo 



27r,5(0i+ 02 +03 + 04)- J2 (ZiZ 2 Z 3 Z 4 ) k (83) 
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M(n 1 ,n 2 ,n 3 ,n 4 ) thus becomes 

oo 

M(m, n 2 , 713,714) = M(n 1 ,n 2 ,n 3 ,n 4 ;k) (84) 



fc=— 00 

where : 

/• 2 ' r # 1 /■ 2w d0 3 f 2 ^ 

M(n 1 ,n 2 ,n 3 ,n 4 ;fc) = — ^ — ^ — ^ — 

(jlvi ■ *T + ^j ■ {Z x Z 2 Z 3 Z 4 ) k . (85) 
The Fourier expansion (9) implies the following integration rule 

(yx m )-ZJ = A(m,j). (86) 



/ 

Jo 



10 2tT 

The calculation of (85) thus becomes straightforward and does not induce any 
summation. The quantity M{n\,n 2 ,n 3 ,n 4 ]k) comes out as a huge expression with 
201 terms, each of them being a product of four hypergeomctric functions. Due 
to the definition (84) where k runs from —00 to 00, and to the symmetry of (84) 
under the permutation of the n^s, these 201 terms reduce to only 16 terms and 
M(ni, n2, 77,3, 714; A;) reads: 

M(ni, n 2 , 713, 714; k) = A(ni,k)A(n 2 ,k)A(n 3 ,k)A(n 4 ,k) 

- A(m, k)A(n 2 , k)A(n 3 ,k)A(n 4 , k + 4) 
A(m, k)A(n 2 ,k)A(n 3 , k + 2)A(n 4 , k + 2) 

+ 2A(m,k)A(n 2 , k)A(n 3 , k + l)A(n 4 , k + 3) 

- 3A(m,k)A(n 2 , k + l)A(n 3 , k + l)A(n 4 , k + 2) 

- A(m, k)A(n 2 , k)A(n 3l k + 3)A(n 4 , k + 5) 
+ A(m, k)A(n 2 , k)A(n 3 ,k + 4)A(n 4 , k + 4) 

- 2A(n 1 ,k)A(n 2 , k + l)A(n 3 , k + 3)A(n 4 , k + 4) 
+ 2A(m,k)A(n 2 , k + l)A(n 3 , k + 2)A(n 4 , k + 5) 
+ 2A(m,k)A(n 2 , k + 2)A(n 3 , k + 3)A(n 4 , k + 3) 

- A(m, k)A(n 2 , k + 2)A(n 3 , k + 2)A(n 4 , k + 4) 
+ A(m, k)A(n 2 , k + 2)A(n 3 , k + 4)A(n 4 , k + 6) 

- A(m, k)A(n 2 , k + 2)A(n 3 , k + 5)A(n 4 , k + 5) 

- A(m, k)A(n 2 , k + 3)A(n 3 , k + 3)A(n 4 , k + 6) 
+ 2A(m, k)A(n 2 , k + 3)A(n 3 , k + 4)A(n 4 , k + 5) 

- A(m, k)A(n 2 , k + 4)A(n 3 , k + 4)A(n 4 , k + 4) (87) 

which is equal to (85) for summation purposes. 

Finally, collecting (81, 84), \^ can be written as : 

OO OO OO OO 

^(4) = i6 W i6 . J2 }2 }2 ^" i+n2+n3+ " 4 x 

m— 712—0 ri3— ri4— 

C(ni,n 2 ,n 3 ,n 4 ) • M(m, n 2 , n 3 , n 4 ) (88) 
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8. Appendix B 

In this Appendix, we give the explicit derivation of the expression (14) given in Section 
2. 

Considering as in (13), we define the following relations 

II f— V = P ^34 - (1 ~ 3) - (2 ~ 3) (89) 

11 T^it = P ^ - (1 - 3) - (2 - 3) (90) 
where P X2 }_ 3 4 is given in (15). Using the constraint (6), P[ 2 l 3i reads : 

A ( 2 Z -34 = + {Z 2 - Z X ){Z Z - Z 4 )(l + ZiZ 2 + Z 3 Z 4 ) (91) 

with = (Z 2 - Z X ){Z 3 - Z 4 ) Zl ^ 2 2 (92) 

(I-Z1Z2) 2 

Note that, while and are completely symmetric under the permutation 

of the angles (pi, the right hand side of (89,90) are anti-symmetric under the same 
transformation. 

From (13, 89, 90) and using the symmetry of the rest of the integrand (namely 
i?( 4 ) times the product over yi variables), we can write as: 

H {A) = \p[ 2 1 m • (^34 - (1 ~ 3) - (2 ~ 3)) (93) 

We obtain given in (14), from (93) with the first term at the right hand 

side of (91). Again using the property of complete symmetry of the integrand, that 
part of (93) with the second term at the right hand side of (91) can be written as: 

1 



(x) 

12-34 



((Zi - Z 2 )(Z 3 - Z 4 )(l + Z X Z 2 + Z 3 Z 4 ) 



- (1 3) - (2 3)) (94) 
and is identically null. 

9. Appendix C 

In this Appendix, we give the explicit integration rules corresponding to the form (16) 
given in the main text. For this purpose, let us consider 

(u^^y^ (95) 

where is given in (92). Expression J(ni, n 2 , 713, n 4 ) has the following properties: 
J(n, n, ri3, n 4 ) — J(ni,n 2 ,n,n) = 

J(ni,n 2 ,n 3 ,n 4 ) = - J(n 2 , m, n 3 , n 4 ) = J(n 3 , n 4 , m, n 2 ) (96) 
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where the last identity comes from the fact that 
Z\Z 2 Z 3 Z 4 



(97) 



(1-Z 1 Z 2 ) 2 " (\-Z 3 Z 4 f 

Changing the integration variables from (<f>i, (f> 2 , (j> 3 ) to (<j> = + (t> 2 , 4> 2 , 4> 3 ) , A^\ 
defined in (92), becomes: 

1 



1(4) 



(l~Zf 



(J— + Z 2 Z 3 )-Z-^ -^z' 
\Z 2 Z 3 J z 3 Z2 



(98) 



with Z = Z X Z 2 . Taking in the form (98), expanding j/ix™ 1 and y 4 x^ 4 , and using 
the Fourier expansion (9), the integration over <p 2 and (j) 3 is straightforward! |. Some 
standard manipulations give: 

r 2iv 11 00 00 



J(n 1 ,n 2 ,n 3 ,n 4 ) = - — ■ E E 

A — — . ^ A — — ^ 



D(m,n 2 ;i) D(n 3 ,n 4 ;j) ■ (99) 

(1 - Z) 

where : 

D(ni,n 2 ;k) = A(m, k + 1) A(n 2 , k) - A(m, k) A(n 2 , k + 1). (100) 
From the definition (10) of A(n, k) , one can easily show that 

£>(ni,n 2 ;-fc) = -D(m,n 2 ;k- 1) (101) 
With this property (101), and after some manipulations, J(ni, n 2 , n 3 , n 4 ) reads 



00 00 



J(ni,n 2 ,n 3 ,n 4 ) = ^2^2D(n 1} n 2 ;i) D(n 3 ,n 4 ;j) N(i,j) (102) 

i=0 j=0 

where 

AJ. 1 _ Z^+l 1 - Z 2j+1 



N ^ = J f« z ^^z-^z- « 

which, by expansion in powers of Z and integration, simply reads: 

N(i,j) = 2min(i,j) + l. (104) 

Finally, taking the form (102), together with (104), and using the definition (10) of 
A(n, k), we obtain: 

00 00 

J( ni ,n 2 , n 3 , n 4 ) = ^ E E ^ 3+4k+2 + 1) (1 + 9 (j - 1)) • 

j=0 k=Q 

(d(ni,n 2 ; k) d(n 3 ,n 4 ;j + k) + d(ni,n 2 ;j + k) d(n 3 ,n 4 ; k)j (105) 
where d(ni,n 2 ; k) is defined in (20). 

|| see the integration rule (86). 
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10. Appendix D 

In this Appendix, we give all the necessary details of the derivation of the result (17) 
given in the text. We consider x^ m the form (2), with given in (14). 

The first step is to consider the expansion of the product ■ Pi 2 -s4- Let us 
give the expansions, in the x^s variables, for the first two terms in the brackets in 
(15): 



l- Xl x 2 l-x 3 x 4 ^ Q ^ o 

9(m -1)9(1-1) (2 min(m, I) -I) 



(x 1 x 2 ) 2 _ _ 

2 

6{m-l- 2) 6>(m - 2) (1 + 9(1 - 1)) . 



i? (4) • = EE (*!*»)*" ( W • (107) 

1 ' Tl ' T2 m=0 Z=0 



' ^12-34 can be expanded as 

R (4) ■ = (X1~X 2 ). (X 3 - X 4 ) ■ E E w2m+2 ' +2 ' ( 108 ) 



oo oo 



m=0 Z=0 

(X1X2)" 1 [x^Xif ■ c(m, 



with 



c(m, Z) = 0(m - 1) 0(Z - 1) • (2 min(m, Z) - 1) 

+ 9(m - I - 2) 0(ra - 2) • (1 + 0(Z - 1)) 
+ 0(Z - m - 2) 9(1 - 2) • (1 + 0(m - 1)) . (109) 
Considering x in the form (2), with H^' given by (14), taking the expansion of 
kl in the form (108), performing the angular integration using the definition 
(95), and collecting all terms, gives 

00 00 

= 16 • E E ^ 2m+2 Pc(m,p)I(m,p) (110) 

m— p— 

where : 

I(m,p) = J(m, m + l,p,p + 1) — J(m,p, m+ l,p+ 1) 

— J(m, p + l,p, m + 1). (HI) 
Using the identities (96) for expression J(ni, n 2 , TI3, 714), one has : 
J(p,m) = I(m,p) 

I(m,m) = I(m,m+ 1) = 0. (112) 

The symmetry of c(m,p), together with the relations (112), allow us, after some 
manipulations, to write x^ as 

00 00 

X (4) = 16 w w • E E ™ 4m+2 " (2™ + 1) • /(m, n + m + 2) (113) 

m— n— 

where we have used the following identity c(m, n + m + 2) = 2m + 1. 
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From the previous equation, taking the definition (111) together with the 
expression (105), one obtains 

oo 

X W = 16w 12 • J2 w 4m+4k+2n+2 ^ (2m + 1) (2k + 1) x 

m,k : n,j— 

\{l + 6(j-l))-V(m,k,n,j-2) (114) 

where V(m,k,n,j) is defined in (19), with the constraints : 

V(m,k,n,-2) = V(m,k,n,-1) = (115) 
Thus, reads : 

x w = i6 w w ■ jr w 4 ™+ 4fe + 2 «+ 2 J x 

m,k : n,j— 

(2m+l)-(2k + l)-V(m,k,n,j). (116) 

The last simplification comes from the following identity on V(m, k,n,j) 

V(k,m,j,n)= V(m,k,n,j). (117) 

which is obtained with the help of the symmetry of a(n, k), namely a(n, k) — a(k, n). 
Finally, using the general identity 



oo oo oo oo 



£ 5>(P1'P2)= (1 + ^-1)) X 

Pi —0 p2— p-0 i-0 

^(s(f> + i,p) + s(p,p + *)) (118) 

and the identity (117), the expression (116) for \^ can, after some manipulation, be 
written in the form (17) displayed in the text. 

11. Appendix E 

P w (x) = 192598769664000- 943722860380160 x 

+ 1154055263764480 x 2 + 223612469147648 x 3 
+ 498965092419008 x 4 - 4709824359388336 x 5 
+ 6098813329440179 a; 6 - 2687506699337617 a; 7 

- 752969324018818 x s + 1919011581320339 a; 9 

- 1526656430056013 a; 10 + 660280621356073 a; 11 

- 134468923815612 a; 12 + 7980003107181 a; 13 
+ 671991155376 a; 14 - 17122807680 a; 15 

- 5111390208 a; 16 + 43929600 a; 17 (119) 



P 9 (x) = 39290149011456000- 297447362118287360a; 

+ 785045342889902080 a; 2 - 738039510467911680 a; 3 
+ 115430038068948224 a; 4 - 1014554971757285568 a; 5 
+ 3530827584228091380 a; 6 - 4191676869770939784 a; 7 
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+ 2038314404276922993 x 8 + 302242744021242129 x 9 

- 1228406406705470056 x 10 + 1043856642401535101 x 11 

- 506114460611418219 x 12 + 139123852417110463x 13 

- 18756086615529738 a; 14 + 585977956065027 x 15 
+ 80289868813296 x 16 + 338076239232 x 17 

- 554882826240 x 18 + 4524748800 x 19 (120) 



P 8 (x) = - 197221140135936000+ 2237384185085952000 x 

- 9353251554131968000 x 2 + 17417661939469352960 x 3 

- 12764493479941255168 x 4 + 5721656503410538176 x 5 

- 31513487740861948240 x 6 + 74327124872294722259 x 7 

- 74793556916527926194 x 8 + 31202576844037093699 x 9 
+ 8393285642528898691 x 10 - 21583664946799604164 x 11 
+ 16831944005767359920 x 12 - 7706219890490209273 x 13 
+ 2036277808796847523 x 14 - 267301496692115881 x 15 

+ 8265260736796220 x 16 + 1095242827496832 x 17 

+ 4127794059264 x 18 - 7320844222464 x 19 

+ 5 85 1 42 2 7 2 00 x 20 (121) 



Prix) = 764231918026752000- 16971026878221516800 x 
+ 115553546628494786560 x 2 - 357128015448284200960 x 3 
+ 527925844930277396480 x 4 - 327328761204753779200 x 5 
+ 264654998297650108784 x 6 - 1129192842855079731378 x 7 
+ 2154212025640162544967 x 8 - 1923633718567453943246 x 9 
+ 708245350538757849333 x 10 + 257335760183269925607 x 11 

- 530877971609916541978 x 12 + 387828430484350338476 x 13 

- 169278351316199852707 x 14 + 43213759197299807177 x 15 

- 5540993509581072863 x 16 + 171042063471920404 x 17 
+ 21596554335396864 x 18 + 67975771296768 x 19 

- 139524297965568 x 20 + 1091386982400 x 21 (122) 



P 6 (x) = 2366653681631232000+ 15657511469396787200 x 

- 418375361108283228160 x 2 + 2413740972333717913600 x 3 

- 6304389540169208659968 x 4 + 7885479774610742167552 x 5 

- 4194762998349204963264 x 6 + 4787870532372577669360 x 7 

- 18347979159424299125323 x 8 + 30953188609834801509491 x 9 

- 25453291295484088365791 x 10 + 8618062554065795648668 x 11 
+ 3501102053489313155231 x 12 - 6497943932928673199052 x 13 
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+ 4524991429929970477637 a; 14 - 1899642804937798765560 a; 15 

+ 470903276066137118684 a; 16 - 59187343449427131707 a; 17 

+ 1841781629480835082 a; 18 + 217686233501214240 a; 19 

+ 505738155562752 x 20 - 1358525988175872 a; 21 

+ 10379070873600 x 22 (123) 



P 5 (x) = -4240254512922624000 +62906863287743283200 a; 

- 176844924247909335040 a; 2 - 926585736920055152640 x 3 

+ 6926618061828909694976 x 4 - 17270267632607743877120 x 5 
+ 19088006038191340885888 x e - 7801708756272263419344 x 7 
+ 10929800840921163293912 x 8 - 45123693732451251946951 x 9 
+ 72469347127122031920503 x 10 - 57004940985415997264335 x 11 
+ 18768946563950160696028 x 12 + 7071447265880540507399 x 13 

- 12955437249853164400164 x 14 + 8753288431922263351485 x 15 

- 3567263867581256107712 x 16 + 863585677517038259840 x 17 

- 106841634606075117511 x 18 + 3388641803935920834 x 19 
+ 366870695477680992 x 20 + 460940531245824 x 21 

- 2210964576190464 x 22 + 16454358835200 x 23 (124) 



P 4 (x) = 2366653681631232000- 121147207913255731200 x 
+ 1204116891510820044800 x 2 - 4471702378853957632000 x 3 
+ 3882247113210548060160 x 4 + 19050981728126452531200 x 5 

- 60173463653707794252800 x 6 + 61647224332324323946560 x 7 

- 8902468613107108100784 x 8 + 12826049043414461768619 x 9 

- 130168614999998943881548x 10 + 224001261613376678757254X 11 

- 177856326190823713976989x 12 + 61054390293590247070213x 13 
+ 16702537875405165273835 x 14 - 34504858377234313267137 x 15 
+ 23145136640956770024901 x 16 - 9263503411249994637068 x 17 
+ 2206094527520924501173 x 18 - 270104745702122623849 x 19 

+ 8842661536198397940 x 20 + 853061197515225600 x 21 

- 104681638038528 x 22 - 4965596931735552 x 23 

+ 35863949721600 x 24 (125) 



P 3 (x) = -11833268408156160000+ 296345375978029056000 x 

- 2144074027349363916800 x 2 + 6754783276113459937280 x 3 

- 8346802888116933918720 x 4 - 3858851854352328468480 x 5 
+ 18632235498956809831680 x 6 - 5896832882414934369872 x 7 

- 20092581953702558202578 x 8 + 9039258514198720200789 x 9 
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+ 30407815233273202124059 a; 10 - 44070504829347073747333 a; 11 
+ 24353838992694268159678 a; 12 - 2634376094447814881673 a; 13 

- 5786419996453628346078 x 14 + 5219893433655665719185a; 15 

- 2420429903132324343712 a; 16 + 639931684854726532004 a; 17 

- 85290635341359290409 a; 18 + 3288947218141266164 a; 19 
+ 249618983188012416 a; 20 - 1896484503705600 x 21 

- 1406458574438400 a; 22 - 98 7 981120 00 x 23 (126) 



P 2 (x) = 36978963775488000000- 561176784419998924800 x 
+ 3011344824723628359680 a; 2 - 7448695810925636485120 a; 3 
+ 7522952469866279936000 a; 4 + 1571640457415428633600 a; 5 

- 6448286163781893470272 a; 6 - 5723224826762698171248 a; 7 
+ 17523982043162827496339 a; 8 - 10871737924438446047496 a; 9 

- 1585720216853837551908 a; 10 + 4808127638071252147363 a; 11 

- 1996433111317397547653 a; 12 - 255305750980018902593 a; 13 
+ 746533058733318361745a; 14 - 451215245481551846087a; 15 
+ 139082951228390470494 a; 16 - 20614123445854835699 a; 17 
+ 931298896466844519 a; 18 + 53593802936851536 a; 19 

- 845315534083200 a; 20 - 297344753664000 a; 21 

+ 2006484480000 a; 22 (127) 



Pi (a;) = -14878398029915750400+ 132065832050507120640 a; 

- 489750894279960821760 a; 2 + 843938312965142528000 a; 3 

- 436158709746751846400 a; 4 - 504050378774640931456 a; 5 
+ 362705915662042967696 x 6 + 930077072222989626372 a; 7 

- 1516628759565419342933 a; 8 + 919766159750634600991 x 9 

- 236614768621801865201 x 10 - 4425339208940491244 a; 11 
+ 15751985529081115811 a; 12 + 356035582683101460 a; 13 

- 3725421245175199951 x 14 + 1704676309420947212 x 15 

- 304427228661108302 x 16 + 16257997137256137x 17 
+ 635239438099728 x 18 - 16870300732800 x 19 

- 3641326080000 x 20 + 23063040000 x 21 (128) 

12. Appendix F 

Let us recall a few elementary facts on factorization of differential operators, without 
introducing any sophisticated P-module theory. Two differential operators T>i and 
2?2 of order TV are said to be equivalent (see for instance equation (5) in [21]) if there 
exist two (intertwining) differential operators R and S of orders at most N — 1, such 
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that : 

V 2 -R = S-Vl (129) 

The factorization of differential operators is not in general unique. There can be 
several factorizations, even an infinite number of factorizations, as can be seen in the 
simple exemple : 

d 2 (d a \ f d 



dx \dx ax + 6 v 
where the factorization (130) holds for every constant a and b. This can be seen as a 
consequence of the relation : 

which means that the multiplication by ax + b transforms the solutions of the 
differential operator d/dx into a solution of d/dx — a/(ax + b). Relation (131) 
is a simple exemple of equivalence (129) of two differential operators. The occurrence 
of this infinite set of factorizations can be seen to be a consequence of the fact that 
the various differential operators of (130), (131) have solutions the ratio of which 
can be rational functions. In general this is not the case and one will only have a 
finite number of factorizations, these various factorizations being basically the same 
up to some permutation of the differential operators, each operator in the factorization 
being changed ("by some transmutation" (129)) into an equivalent one during the 
permutation process. 

This is typically the situation we encounter with the operator £ w associated with 
the ordinary differential equation of order ten satisfied by ■ 

This order ten operator Ciq has 36 factorizations with two operators of order four 
(Mi and M 2 ), ten operators of order two (No, N\, ■ ■ ■ , iV 9 ) and 26 operators of order 
one (Lq, Li, • • • £25)- All these operators have been normalized in such a way that the 
coefficient of the highest derivative is normalized to +1. The operators Lq, Li, L2, 
N and N$ below are such that L (Sq) = 0, Li(<Si) = 0, L 2 (S 2 ) = 0, N (S 3 ) = and 
N^(S*) = 0. These 36 factorizations* of C w read 

-ClO = -^8 ' M2 ' L25 ■ L\2 ■ L3 ■ Lq 

= Mi • 7V 9 . L 25 ■ L 12 ■ L 3 ■ L a 

= Mi • L24 • JV4 ' L12 ■ L3 ■ Lq 
= Ml • L24 ' L13 ■ Nq ■ L3 ■ Lq 
= Ml • Zv24 ' Ll3 ' Ln ■ N3 ■ Lq 

= Mi • L 2 ± ■ L 13 ■ L 17 ■ Ln ■ Nq (132) 



C 



10 



N 8 - 


M 2 • 


L25 ' 


L14 ■ 


L 4 - 


Lq 


Mi 


■Ng- 


L 2 5 ' 


L\4 ■ 


L^ 


Lq 


Mi 


■ L 2 4 


• N A ■ 


L14 ■ 


L 4 - 


Lq 


Mi 


■ L 2 4 


■ Li 5 


■ N 7 ■ 


L 4 - 


Lq 


Mi 


■ L 2 4 


■ Li 5 


■ L w 


■N 3 


■Lq 


Mi 


■ L 2 4 


■ Li 5 




■L u 


■Nq 



(133) 



* The DFactor command [27] of DEtools can give the second factorization in (132), namely £10 
M\ ■ Ng ■ L25 • L12 • L3 ■ Lq. 
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£10 — Ng • M2 ■ L25 ■ Lis ■ L$ ■ L\ 
= M 1 -N 9 - L 25 ■ L m ■ L 5 - L 1 
= Mi • L 24: ■ N 4 ■ L 1S ■ L 5 • L 1 
= Mi • L 2 4 ■ L 19 ■ N 5 ■ L 5 ■ L x 
= Mi ■ L2A ■ L\g ■ L23 ■ N± ■ L\ 

= Mi • L 2i ■ L 19 ■ L 23 ■ L w ■ N (134) 



N 8 - 


M 2 • 


L25 • 


L14 ■ 


Lg ■ Li 


Mi 


■N g - 


L25 • 


L14 ■ 


Lq ■ L\ 


Mi 


■ L 2 4 


■N 4 - 


L14 ■ 


La ■ L\ 


Mi 


■ L 2 4 


■ L\s 


■N 7 . 


Lq ■ L\ 


Mi 


■ L 2 4 


• Li 5 


• L 2 o 


■ N\- L\ 


Mi 


■ L24 


■ L15 


■ L 2 o 


■ L w ■ No 



(135) 



N 8 - 


M 2 • 


L25 ' 


Lm ' 


L 7 - 


L 2 


Mi 


■N 9 - 


L25 • 


Lis • 


L 7 - 


L 2 


Mi 


■ L 2 4 


• N A ■ 


Lis • 


L 7 - 


L 2 


Mi 


■ L24 


■ Lig 


■N 5 - 


L 7 - 


L 2 


Mi 


■ L24 


■ Lig 


■ L21 


■N 2 


■L 2 


Mi 


■ L 2 4 


■ Lig 


■ L 2 i 


■L 9 


■N 



(136) 



N 8 - 


M 2 • 


L25 ■ 


L12 ■ 


Ls- 


L 2 


Mi 


■Ng- 


L25 ' 


L12 ■ 


Ls- 


L 2 


Mi 


■ L 2 4 


• N A ■ 


L12 ■ 


Ls- 


L 2 


Mi 


' -^24 


' -^13 


■No- 


Ls- 


L 2 


Mi 


' ^24 


' £l3 


■ L22 


■N 2 


•Li 


Mi 


• -^24 


• £l3 


■ L 22 


■L 9 


■N 



All these 36 factorizations can be seen as the consequence of "transmutations" like 
(129), where one differential operator becomes an equivalent one, thus yielding some 
elementary permutation of the operator. For instance, the two differential operators 
of order four, namely Mi and M 2 , are equivalent, their "intertwiners" being two order 
two differential operators N s and N 9 : 

N 8 ■ M 2 = Mi • N 9 (138) 

Actually, one sees clearly that the two first factorizations in (132) are just deduced 
from one another by (138). 
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Nine order- two differential operators N ,---,N S can be seen to be related by 
various equivalence relations requiring the introduction of 26 differential operators 
LiS of order one (the intertwiners) : 



(139) 



Nl 


Li 


— L10 


N , 


N 2 


L 2 = 


L 9 -No 


N 3 


L 


= L n 


N , 


N 4 


L\2 = 


-- L 13 ■ N 6 


N 4 


Lu 


= L15 


N 7 , 


N 4 


Lis = 


= L 19 ■ N 5 


N 5 


L 5 


= L23 


N lt 


N 5 


L 7 = 


L 2 i ■ N 2 


N 6 


L 3 


= L 17 


N 3 , 


N 6 


Ls = 


L22 • N2 


N 7 


L A 


= Liq 


N 3 , 


N 7 


L 6 = 


L 20 • Nx 


N 9 


L25 


= L 2 4 


N 4 . 









The 26 order-one differential operators L^s can also be seen to be related by 
various equivalence relations. The L^s being order one differential operators, their 
intertwiners (R and S in (129)) should be order zero differential operators, that is 
some function fij. As a consequence of the normalization to +1 of the highest order 
derivative (here d/dx), the two intertwiners R and S in (129) are identical and one 
thus gets relations like : 

Li ■ f itj - fij ■ I.,, (140) 

Furthermore one also verifies the following relations between the order one 
operators Li's : 

(141) 



(142) 



(143) 

It is straightforward to see that a relation like Lj • Lj = Lk ■ L m yields relations like : 
Li ■ (Lj — L rn ) = (Lk — Li) ■ L m , 
(Li — Lk) ■ Lj = Lk ■ (L m — Lj), 

fi.m = Lj — L rn = Lk — L; t , fkj = L m — Lj = Li — Lk 

The number of intertwiners f±j between the operators of order one in (140), is at first 
sight 325, but due to the above identities only 15 are independent. 
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